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The surface tension a of the confined-deconfined interface is calculated in 
pure SU(3) lattice gauge theory at finite temperatures employing the oper- 
ator and integral methods on a lattice of a size 8 2 x N z x 2 with N z = 16 
and 40. Analyses of non-perturbative corrections in asymmetry response 
functions strongly indicate that the use of one-loop values for the response 
functions lead to an overestimate of a in the operator method. The operator 
method also suffers more from finite-size effects due to a finite thickness of 
the interface, leading us to conclude that the integral method yields more 
reliable values for a. Our result with the integral method cr/T c 3 = 0.134(16) 
is consistent with earlier results and also with that obtained with a transfer 
matrix method. Result is also reported on a obtained on a lattice 18 2 x 48 x 4 
with the integral method. 



1 Introduction 



Numerical simulation of pure SU (3) gauge system on a lattice has shown that 
the system undergoes a first order phase transition from a confined phase at 
low temperatures to a deconfined phase at high temperatures 0. The two 
phases can therefore coexist at the transition temperature T c , separated by 
an interface. 

A basic parameter characterizing the interface is the interface tension a. 
A number of numerical work has recently been carried out to determine its 
value, mainly for a system with the temporal lattice size N t = 2. Kajantie, 
Karkkainen and Rummukainen|2|, |3| developed an operator method for mea- 
suring the tension and reported the value cx/T c 3 = 0.24(6) for N t = 2[||]. 
Independently, Potvin and RebbiQ proposed to employ an integral of the 
derivative of the free energy in the parameter space of the coupling constant 
to measure the tension, and found cr/T c 3 = 0.115(13) for the same temporal 
lattice size|J. A factor two discrepancy between these results has motivated 
further studies of the interface || ^, |[ [|, which yielded values of cr/T c 3 similar 
to that of the integral method. The method of histograms [ID| and the tech- 
nique of transfer matrix|| used in these studies are quite different, however, 
from the operator and the integral method. Hence the discrepancy of the 
original results obtained with the two methods has not been resolved. 

In this article we report on our attempt toward an understanding of the 
origin of the discrepancy through a comparative study of the two methods. 
Since a possible cause of the discrepancy is the use of a different spatial 
lattice size (8 2 x 40 x 2 for the operator method || and 8 2 x 16 x 2 for the 
integral method|| ), we carry out simulations for both sizes of lattice with 
both methods. We also increase statistics by a factor of two to four over 
those of ref . |3], [| . We find that a part of the discrepancy can be understood 
by the effect of the finite thickness of the interface. Comparing the results on 
a large lattice with N z = 40, where this problem is absent, we find, however 
that the discrepancy is not completely removed for N t = 2. 

A nice feature of the integral method is that it directly relates plaquette 
expectation values to the interface tension and no other input is needed. 
In contrast the operator method requires the values of response functions 
of the coupling constant under an anisotropic deformation of lattice. Up 
to now they have been estimated only in one-loop perturbation theory O], 
whose applicability may be questioned for small values of the inverse gauge 
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coupling constant (3 = 6/g 2 , such as the critical value (3 C for N t = 2. We 
discuss the effect of nonperturbative corrections to these coefficients. 

Another problem which became apparent in recent work on the interface 
tension is the difficulty of measuring its value for lattices with the temporal 
size N t > 4. For the case of N t = 4 an initial attempt with the integral 
method failed to obtain a non- vanishing value 0. An alternative integral 
method employing an external field coupled to the Polyakov loop was then 
devised, and a value a/T c = 0.027(4) has been obtained^]. We attempt a 



high statistics measurement of the interface tension with the original integral 
method on a 18 2 x 48 x 4 lattice. The results will also be reported in this 
article. 

This paper is organized as follows. In Sec. 2 we briefly recall the operator 
and the integral method and describe our simulations. Our results for the 
interface tension are presented in Sects. |3] and |] for N t = 2 and 4, respectively. 
Sec. H is devoted to conclusions. 



2 Method and simulation 

We consider the pure SU (3) gauge system with the standard plaquette action 
on a lattice of a size N x x N y x N z x N t with the size N z in the z-direction 
at least twice as large as those in the x- and y-directions. If one splits the 
lattice into two halves in the z-direction and chooses (3 in one half of the 
lattice to be slightly below the critical value /?i = /3_ = f3 c — 5/3 and that 
of the other half slightly above @ 2 — P+ = Pc + 5/3, the system in the first 
half will be in the confined phase and the other half in the deconfined phase. 
With the periodic boundary condition employed for the gauge field in all 
four directions, there will be two interfaces in the x — y plane, each with a 
transverse area A = N x N y a 2 with a the lattice spacing. 

The operator method^], [| starts from the following thermodynamic re- 
lation defining the interface tension a in terms of the free energy F, 



ldF 
a= 2dA 



(1) 



T,V 

where the factor 1/2 is to take into account the presence of two interfaces. 
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For the standard plaquette action, working out the derivative yields ||, 
a A 1 

~7fT = 7y *^2lPn + 3(Ct(f3 n ) — C s (/? n )](2i 3 „ ; i2— P n; 13 — -Pn;23 — 2P n; 03 + -P ra; 01 + -Pn;02) 
* n 

(2) 

where T = 1/ (iVja) is the temperature, fa is the local value of (3 at site n 
taking either of the values fa = fa or fa = P+, and P n]fU/ = (|ReTr[/ n;/LU/ ) is 
the plaquette expectation value in the (//, z/) plane at site n. The functions 
c s (/3) and c t {(3) are the response functions of the gauge coupling constant with 
respect to an anisotropic deformation of the lattice in the spatial and tem- 
poral directions: q = (dg^ 2 /d^)^ = %, i = s, t, and £ = a s /a t . In perturbation 



theory these functions are expanded as|TT 



c s {(3) = 0.20161 + 0(/T 1 ) , 

dp) = -0.13195 + Oip- 1 ) . (3) 
The integral method [5|] is based on the relation 

°A = ~[F(^,f3 + )-~{F((3 + ,f3 + ))+F(fa,fa)}\ 

= ~[{F((3_,(3 + )-F(fa,{3„)}-{F(fa,(3 + )-F(fa,P + )}\ (4) 

where F(fa, fa) is the free energy of the system with the two halves having a 
coupling fa and fa. Denoting expectation values for such a system as (•) l a 1( g 2 
and decomposing the total action into the contributions of the two halves 
S = faS\ + faSz, we have 

l dF(fa,fa) 

j, -qq = -{Si/Pifh » i = l,2 (5) 

Expressing the free energy differences in (|J) through integrals of (|5]) we obtain 
a A 1 rP+ 

-T = l [ -J p tbWfi-fh+h dp 1 (S 1 )^ + ). (6) 

For both the integral and operator methods numerical simulations yield 
an estimate of a for a finite difference 5f3 between the couplings fa = fa — 5(3 
and fa = fa + 5/3 chosen for the two halves of the lattice. The physical value 
of the interface tension is obtained by an extrapolation to the limit Sf3 — > 0. 
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Our simulations are carried out on a lattice of a size 8 2 x N z x 2 with 
N z = 16 and 40, and also on a 18 2 x 48 x 4 lattice. Gauge configurations are 
generated through the standard pseudo-heat bath algorithm with eight hits 
per link. 

In order to create an interface we take the gauge coupling constant to be 
/?i for all the plaquettes starting at the site (x,y,z,t) with < z < N z /2 
and pointing in the positive \xv direction. For the rest of plaquettes at 
N z /2 < z < N z , we assign fi 2 - Thus the coupling constant jumps sharply at 
z = and N z /2. This assignment is slightly different from that adopted in 
ref. |3|, [| where the lattice is split into two at z = 1/2 and N z /2 + 1/2 and 
the coupling (3\ or /3 2 is used for updating links in the left or right half of 
the lattice, while the mean value (/3i + ^2)j2 is chosen for links connecting 
z = to 1 and N z /2 to N z /2 + 1. We prefer our assignment since (i) thermo- 
dynamic interpretation is clearer with (3 associated to plaquettes so that the 
convergence to thermal equilibrium is guaranteed in numerical simulations, 
and (ii) there is no ambiguity in expressing dF/df3i (i = 1,2), needed in 
the integral method, in terms of plaquette expectation values. We note that 
the physical value of the interface tension obtained in the limit /?i — 02 — > 
should not depend on the details of assignment. This is found to be the case 
as will be shown below. 

Another choice one has to make in applying the two methods is the value 
of the critical coupling (3 C) which involves ambiguities on a finite lattice. 
A possible definition of (3 C on a finite lattice is the peak position of the 
susceptibility of the Z(3)-rotated Polyakov loop. Carrying out runs on a 
8 2 x N z x 2 lattice with f3i = /3 2 and applying the spectral density method |i~3f 
to locate the peak position, we find that (3 C — 5.09146(63) and 5.09409(38) for 
N z = 16 and 40. These values are slightly higher than the choice j3 c = 5.09 
used in refs. @, 3- In order to facilitate a comparison of results, however, 
we have decided to use the same value (3 C = 5.09 for both N z = 16 and 40 
lattices with N t = 2. For N t — 4 we use the estimate (3 C = 5.6924 obtained 
on a 24 2 x 36 x 4 lattice PR. 



In Table |TJ - |2| we list the parameters of our runs on N t = 2 lattices 
together with the plaquette expectation values obtained. In Fig. |] the cor- 
responding points in the f}\ — j3 2 plane are plotted. In these tables the first 
group of runs with (3\ = 02 corresponds to the simulation points on the line 
AF of Fig. |l|. The second group of runs located on the line CD are used for 
measuring the interface tension with the operator method. For the integral 
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method we follow ref. |J and modify the integration path as shown in Fig. [I], 
namely the path ABCD is used for the first integration in (||) and DCEF 
for the second one. The results along the paths BC and CE are collected 
in the third and fourth group in Tables p] and |2|. With these runs we can 
measure the interface tension for four values of 5(3 = 0.1, 0.15, 0.2, and 0.3. 
The runs on an N t = 4 lattice are summarized in Table |3|. They cover the 
cases 5(3 = 0.02, 0.03, 0.04, and 0.06. The integration contour is similar to 
those of Fig. |l| The runs for N t = 2 were made on the 16 processor version 



of QCDPAX and those for N t = 4 on QCDPAX itselfp 



For error estimation of the interface tension in the operator method we 
apply the jackknife procedure. An analysis of the bin size dependence showed 
that the estimated error is quite stable over a wide range of bin size of O(100) 
- O(10 000) sweeps. The errors we quote are for the bin size of about 1000 
sweeps. Near the critical point, larger bin size is chosen. Errors of plaquette 
expectation values in Table [3] are also calculated with the same bin size. 

In the integral method we use the first-order formula to numerically in- 
tegrate the plaquette averages. Jackknife errors of plaquette averages are 
combined quadratically to obtain the statistical error of the integral. We 
should note that a finite integration step size causes a systematic error in ad- 
dition to the statistical one. We estimate this error by comparing the results 
of the first-order integration with those of the natural spline fit of the inte- 
grand. In the worst case, the systematic error estimated in this manner is of 
comparable magnitude with the statistical error computed for the first-order 
formula; in most cases, however, the systematic error is much smaller than 
the statistical error. For the final values of errors, we add the systematic 
error to the statistical one. 



3 Results for N t = 2 
3.1 Operator method 

In Table |] we list the results of the operator method for cx/T c 3 for N t = 2 
lattices. The one-loop values (^|) are used for the response functions. Our 
estimate of cr/T c 3 at 5(3 = obtained by a linear extrapolation in 5(3 is also 
given. 

In Fig. |2| we compare our data on an 8 2 x 40 x 2 lattice with those of 
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Kajantie et aZ.[E| on the same size of lattice. We observe that our results at 
a finite 5(3 are larger than their values. The discrepancy is ascribed to the 
difference in the assignment of the coupling constant for creating the interface 
as was discussed in Sect. @. The discrepancy diminishes for a smaller 5j3, and 
our extrapolated value cr/T c 3 = 0.199(34) is consistent with the result 0.24(6) 
by Kajantie et al. 

Our results on an 8 2 x 16 x 2 lattice is smaller than those on an 8 2 x 40 x 2 
lattice as shown in Fig |[ The difference increases with decreasing 5(3, and 
a linear extrapolation toward 5(3 = results in a factor two difference in the 
values of cr/T c 3 for the two lattice sizes (see the last column of Table |]). 

To understand the origin of this size dependence, we study the spatial 
thickness of the interface. In Fig. || we show the z-dependence of the Polyakov 
loop rotated to the nearest Z(3) axis fl TOt (N z ) for various 5(3. For the larger 
lattice with N z = 40 we observe a clear plateau away from the the interface 
centers at z = and N z /2 for all values of 5(3. In contrast, a plateau is barely 
visible for the N z = 16 lattice even for 5(3 = 0.3. Furthermore, for 5(3 ~ 0.15, 
f2 rot (16) deviates not only from the asymptotic values of f2 ro t(40) but also 
from fi rot (40) at the same distance from the interface center even at the 
middle points z = 4 and 12 between interfaces. These observations suggest 
that the half width of the interface region is larger than 4 and consequently 
the two interfaces are affecting each other on the N z = 16 lattice. 

Similar behavior is also found for the ^-dependence of the interface free 
energy, albeit with larger errors. In Fig. [|we plot the quantity a(d) obtained 
by restricting the summation in (^j for the interface tension to the two regions 
within a distance d from the interface centers at z = and N z /2. For 
5(3 = 0.1, we find that a(d) on the N z = 40 lattice saturates only after 
d = 5 — 6, and that a(d) for N z = 16 is clearly not saturated within the 
distance d < 4 available on this lattice. Furthermore the value of a(d) 
for N z = 16 is smaller than that for N z = 40 at d ~ 2 - 4. From these 
observations, we conclude that N z = 16 is not large enough to have two 
free interfaces within the spatial extent of the lattice. We also note that 
the smaller interface free energy observed for N z = 16 corresponds to an 
attractive interaction between opposite interfaces. 

Another problem with the results in Table |] is that they were obtained 
employing one-loop perturbative results for the response functions c t and c s 
and the lattice bare coupling /3 in (0) . This represents a questionable point 
since the critical coupling constant on an N t = 2 lattice is large. In order to 
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examine how non-perturbative effects possibly affect the values of response 



functions [16], we recall that the combination q + c s is proportional to the 



scaling beta function: c t + c s = — | (f^). A non-perturbative estimate of this 



quantity is available from Monte-Carlo renormalization group studies fl7 
which shows a large deviation of the f3 function from the perturbative result 
for (3 ~ 6. The coefficients c t and c s also appear in other thermodynamic 
quantities, the energy density e and the pressure p[(0|. In fact the combi- 
nation e — 3p is proportional to c t + c s , while e + p depends both on c t + c s 
and Ct — c s . The problem with perturbative estimates of response functions 
manifests here not only in a scaling violation of e and p but also in a non- 



vanishing pressure gap at the first order deconfining transition [[ijj. We can 
therefore determine c t — c s at (3 C by combining the results for the beta func- 
tion from Monte-Carlo renormalization group studies with the requirement 
of a vanishing pressure gap. 

Applying this procedure to the data obtained on 24 2 x 36 x 4 and 36 2 x 
48 x 6 lattices @, we find that 

c t - Cs = -0.52(12) (3 = 5.6925 

= -0.37(16) (3 = 5.8936 (7) 

{c t = -0.24(6) and -0.16(8) and c s = 0.28(6) and 0.21(8) for the two values 
of /?, respectively.) These values should be compared with the perturbative 
value, 

c t -c s = -0.3336 + 0(/T 1 ). (8) 

The non-perturbative estimates ([7]) imply that the value for the interface 
tension will be reduced by a factor of 0.88(8) for N t = 4 and 0.98(10) fovN t = 
6 from that with the perturbative value using (|8|). This leads us to expect a 
reduction of the interface tension by perhaps an even larger magnitude for 
N t = 2 if non-perturbatively determined response functions are employed. A 
quantitative estimation of the magnitude of reduction, however, requires the 
value of the effective beta function at the critical coupling (3 C = 5.09 which 
is not available at present. 

We should note that the discussion and the conclusion from Fig. |5| con- 
cerning finite-size effects are not affected by the correction to q and c s be- 
cause the main effect of it is a shift of the overall coefficient of o. 
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3.2 Integral method 

In Table [| we summarize our results for cx/T 3 for the integral method. As is 
seen in Fig. |] the data lie on a straight line with respect to 5(3. The results 
of a linear extrapolation are also given in Table [5[ In Figs. [7] and |8| we plot 
the plaquette expectation values in the two halves of lattice as a function of 
(3 along the integration path for the case of 5/3 = 0.1. The area surrounded 
by the two lines gives an estimate of a. 

As in the case of the operator method we find that our values for N z = 16 
has a different slope from the corresponding results of Huang et al. at finite 
8(3 (compare filled circles and triangles in Fig. P), while the value extrapo- 
lated to 5(3 = is perfectly consistent, namely we find cr/T c 3 = 0.113(6) as 
compared to cr/T 3 = 0.115(13) in ref. ||. 

The results on an 8 2 x 40 x 2 lattice are larger by about 5-10% from 
those on an 8 2 x 16 x 2. The difference increases to 20% after extrapolation 
to 5j3 = 0, which, however, is much smaller compared to the case of the 
operator method where a factor two discrepancy is observed. We also note 
that the size effect is not sensitive to the value of 5/3 (Fig. ^), in contrast to 
the case of the operator method discussed in the previous subsection. 

The origin of this difference of lattice size dependence may be understood 
in the following way. We expect the size effect on the interface to increase 
when both (3\ and j3 2 become close to (3 C while keeping (3\ < (3 C < /3 2 . With 
our choice of the integration path as shown in Fig. [I], the size effect is therefore 
largest at the point C, diminishing toward the points D, B and E. Near the 
points B and E and along the paths BA and EF we expect size effects to 
be small because of the absence of the interface. We now recall that the 
interface tension in the integral method is obtained from the integrals along 
the paths ABCD and DCEF, while that with the operator method is taken 
at a point along the path DC. We then expect the magnitude of size effects to 
be smaller in the integral method since the large size effect near the point C, 
directly affecting results for the operator method, is diluted by contribution 
from other parts of the paths suffering less from the effects. Furthermore 
the fact that the paths BC and CE are common to all of our choices of 8(3 
explains a similar magnitude of size effects seen in Fig. [| 

Comparing the results for N z = 40, we find that the value of cx/T 3 ob- 
tained by the operator method (with the perturbative coefficients) is larger 
than that of the integral method by a factor of about 1.5. We ascribe this 
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discrepancy to the presence of non-perturbative corrections to the response 
functions based on the consideration discussed in the previous subsection. 

We also note that our result cr/T 3 = 0.134(16) with the integral method is 
consistent with a value cx/T 3 = 0.139(4) from a determination with a transfer 
matrix method on 4 4 x 64 x 2 - 8 2 x 128 x 2 lattices 0, while a determination 
with a histogram method on 8 2 x 30 x 2 - 14 2 x 42 x 2 lattices gives a 
slightly smaller value 0.092(4). 

4 Results for N t = 4 

Our comparative analysis of the operator and integral methods for N t = 
2 lattices shows that the integral method suffers much less from various 
systematic errors. We therefore adopt the integral method for our study on 
an N t = 4 lattice. 

A difficulty in an extraction of the interface tension for Nt = 4 is the 
smallness of the value of aa 3 . In fact assuming scaling the result cr/T c 3 ~ 0.1 
for Nt = 2 leads to era 3 ~ 0.0015 for Nt = 4. The actual value appears 
even smaller; the value cx/T 3 = 0.027(4) obtained with the integral 
method using an external field coupled to the Polyakov loop translates to 
aa 3 = 0.00042(6). Thus very accurate results of plaquette expectation values 
for 5(3 much smaller than those employed for the case of N t = 2 will be 
needed. Furthermore it is known that a spatial size of at least 4-5 times 
larger than the temporal size is needed for clear first-order signals on an 
N t = 4 lattice [p0| , fLj}] . This requirement was not met in the first attempt || 
toward an interface tension measurement on an N t = 4 lattice, which failed 
to find a finite value of a employing 8 2 x 16 x 4 and 12 2 x 24 x 4 lattices. 
Based on those considerations we chose to work on an 18 2 x 48 x 4 lattice 
and carried out simulations for 5(3 = 0.02 - 0.06. 

As discussed in the previous section, a large N z is required to remove the 
effect of interface thickness. In Fig. we plot Q rot for 5(3 = 0.02 and 0.04 
as a function of z. We observe a plateau for both cases, albeit somewhat 
limited in range for 5(3 = 0.02. From our experience with the N t = 2 lattices 
we consider that the effect of interface thickness is sufficiently suppressed for 
5/3^0.02 

Our result for a is summarized in Table |6| and plotted in Fig. [K]. Taking 
the three points with 5(3 < 0.04 and extrapolating linearly toward 5(3 = 
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we find cr/T c 3 = 0.0171(113). The central value is smaller than the previous 
estimate a/T 3 = 0.027(4) obtained on a 16 2 x 32 x 4 lattice |TJ. We note 



that a recent study with the histogram method using the data of [TJj on 
lattices 12 2 x 24 x 4 and 24 2 x 36 x 4 also gives a larger value 0.0292(22) 
[ f2~l|j . The error in our result, however, is too large to see if there is a real 
discrepancy. 



5 Conclusions 

We have studied the confinement-deconfinement interface tension in the pure 
SU(3) gauge theory on an 8 2 x N z x 2 lattice with N z = 16 and 40 by ap- 
plying the operator and integral methods. Our results confirm the values 
reported previously when compared at the same parameters. Our study of 
the N z dependence shows, however, that the thickness of the interface causes 
misleadingly smaller values for the interface tension a when estimated on a 
lattice of a small length such as N z = 16. This effect is severer for the oper- 
ator method, which explains part of the discrepancy in the interface tension 
obtained with the two methods. We made some quantitative estimates on 
non-perturbative corrections to the responce function q — c s needed in the 
operator method, and argued that incorporating them would remove the rest 
of the discrepancy. We conclude that the integral method yields more reliable 
values for the interface tension, suffering less from various systematic errors. 
Our result cr/T c 3 = 0.134(16) with the integral method on a large lattice is 
consistent with a determinations with a transfer matrix method ||. 

We have also attempted a measurement of the interface tension on an 
N t = 4 lattice with the integral method. The result is smaller than those 



reported in the literature JT2], It has quite a large error, however, in 

spite of the substantial statistics accumulated with 14 days of fully running 
QCDPAX. The problem originates from a rapid decrease of the discontinuity 
in plaquette averages with an increasing N t . A judicious choice of observables 
having a larger jump across the deconfining transition as was done by Brower 
et al. ||12j| will be indispensable for a determination of the interface tension 
closer to the limit of continuum space-time. 



10 



Acknowledgements 

We are grateful to A. Ukawa for valuable discussions and suggestions and 
continuous encouragement. We also thank the members of the QCDPAX 
collaboration for discussions and their help. This project is in part sup- 
ported by the Grant-in-Aid of Ministry of Education, Science and Culture 
(No.62060001 and No.02402003) and by the University of Tsukuba Research 
Projects. 

References 

[I] For a review of evidence, see, A. Ukawa, in Proceedings of Lattice'89, 
Nucl. Phys. B(Proc. Suppl.) 

[2] K. Kajantie and L. Karkkainen, Phys. Lett. B214 (1988) 595. 

[3] K. Kajantie, L. Karkkainen and K. Rummukainen, Nucl. Phys. B 333 
(1990) 100 

[4] J. Potvin and C. Rebbi, Phys. Rev. Lett. 62 (1989) 3062. 

[5] S. Huang, J. Potvin, C. Rebbi and S. Sanielevici, Phys. Rev. D 42 (1990) 
2864; errata 43 (1991) 2056. 

[6] B. Grossmann, M.L. Laursen, T. Trappenberg and U.-J. Wiese, Nucl. 
Phys. B396 (1993) 584. 

[7] W. Janke, B.A. Berg and M. Katoot, Nucl. Phys. B382 (1992) 649; 

[8] B. Grossmann, M.L. Laursen, T. Trappenberg and U.-J. Wiese, Phys. 
Lett. B293 (1992) 175. 

[9] B. Grossmann, M.L. Laursen, Jiilich preprint HLRZ-93-7, revised, to 
appear in Nucl. Phys. B. 

[10] K. Binder, Z. Phys. B43 (1981) 119; Phys. Rev. A25 (1982) 1699. 

[II] F. Karsch, Nucl. Phys. B205[FS5] (1982) 285. 



11 



[12] R. Brower, S. Huang, J. Potvin and C. Rebbi, Phys. Rev. D 46 (1992) 
2703. 

[13] I.R. McDonald and K. Singer, Discuss. Faraday Soc. 43 (1967) 40; A.M. 
Ferrenberg and R. Swendsen, Phys. Rev. Lett. 61 (1988) 2058; 63 (1989) 
1195. 

[14] Y. Iwasaki, K. Kanaya, T. Yoshie, T. Hoshino, T. Shirakawa, Y. Oy- 
anagi, S. Ichii and T. Kawai, Phys. Rev. D 46 (1992) 4657. 

[15] Y. Iwasaki et al., Computer Physics Communications 49 (1988) 449; Y. 
Iwasaki et al, Nucl. Phys. B (Proc. Suppl.) 17 (1990) 259; T. Hoshino, 
PAX Computer, High-Speed Parallel Processing and Scientific computing 
(Addison- Wesley, New York, 1989). 

[16] For an early attempt for a non-perturbative study, see, G. Burger et al, 
Nucl. Phys. B304 (1988) 587. 

[17] A. Hasenfratz et al, Phys. Lett. B140 (1984) 76; K.C. Bowler et al, 
Nucl. Phys. B257[FS14] (1985) 155; A.D. Kennedy et al, Phys. Lett. 
B155 (1985) 414; R. Gupta et al, Phys. Lett. B211 (1988) 132; J. Hoek, 
Nucl. Phys. B339 (1990) 732. 

[18] J. Engels et al, Nucl. Phys. B205[FS5] (1982) 545. 

[19] J. Engels et al, Phys. Lett. B252 (1990) 625. 

[20] M. Fukugita, M. Okawa and A. Ukawa, Nucl. Phys. B337 (1990) 181. 

[21] Y. Iwasaki, K. Kanaya, L. Karkkainen, K. Rummukainen and T. Yoshie, 
preprint UTHEP-252 (1993). 



12 







run parameters 


Pi 






P2 




p 2 -p 




tot. sw. 


tneim. 
















4. 


79 


01 nnn 

Z1UUU 


i nnn 

1UUU 


o 


373370 


(4:5) 








4 


89 


Z1UUU 


i nnn 
1UUU 


o 


386170 


(37) 








4 


94 


1 uuu 


i nnn 

1UUU 


o 


392992 










4 


99 


OUUUU 


i nnn 

1UUU 





400056 


(46) 








5 


00 


OUUUU 


i nnn 

1UUU 


o 


401590 


(48) 








5 


04 


Qfinnn 

OUUUU 


1 nnn 

1UUU 





407723 


(55) 








5 


07 


30000 


i nnn 





41261 


(10) 








5. 


09 


i nsnnn 
ryouuu 


snnn 

oUUU 





427? 


I 


(15) 








5 


11 


qnnnn 

OUUUU 


i nnn 

1UUU 





45551 


(22) 








5 


14 


^nnnn 

OUUUU 


1 nnn 

1UUU 





46652 


(12) 








5 


19 


^nnnn 

OUUUU 


i nnn 

1UUU 





480807 


(95) 








5. 


24 


qnnnn 

OUUUU 


i nnn 

1UUU 





492618 


(75) 








5. 


29 


91 nnn 

Z1UUU 


i nnn 

1UUU 





503312 


(75) 








5. 


39 


24000 


4000 





521? 


125 


(68) 








4.99 


5.19 


onnnnn 
zuuuuu 


1 nnnn 

1UUUU 


0.407549 


(86) 


0.470238 


'92) 


0.062689 ( 


79) 


4.94 


5.24 


9nnnnn 
zuuuuu 


1 nnnn 

1UUUU 


0.399601 


(45) 


0.482138 




0.082537 ( 


52) 


4.89 


5.29 


1 o^nnn 

1ZOUUU 


1 nnnn 

1UUUU 


0.392557 


(42) 


0.492648 




0.100091 ( 


56) 


4.79 


5.39 


onnnnn 
zuuuuu 


1 nnnn 

1UUUU 


0.379668 


(27) 


0.511044 


'35) 


0.131376 ( 


37) 


4.99 


5.04 


30000 


3000 


0.400363 


(64) 


0.407368 


v60) 






4.99 


5.07 


60000 


3000 


0.400476 


(50) 


0.411974 


J 58) 






4.99 


5.09 


60000 


3000 


0.400610 


(45) 


0.415267 


'67) 






4.99 


5.11 


60000 


3000 


0.400846 


(61) 


0.41970 


'31) 






4.99 


5.12 


60000 


3000 


0.40124 


(10) 


0.42504 


J 77) 






4.99 


5.14 


60000 


3000 


0.40338 


(15) 


0.4459 


; io) 






4.99 


5.16 


60000 


3000 


0.40539 


(14) 


0.45930 


;35) 






5.04 


5.19 


60000 


3000 


0.42676 


(32) 


0.47455 


;i5) 






5.07 


5.19 


60000 


3000 


0.44340 


(31) 


0.47709 


;n) 






5.09 


5.19 


60000 


3000 


0.4519 


(16) 


0.477845 


J 97) 






5.11 


5.19 


60000 


3000 


0.45917 


( 12 ) 


0.478778 


'82) 






5.14 


5.19 


30000 


3000 


0.46769 


(14) 


0.47928 


;n) 







Table 1: Run parameters and plaquette expectation values in each domain 
on the 8 2 x 16 x 2 lattice. These simulation points are plotted in Fig. 1. 
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run parameters 


Pi 






P2 




p 2 -p 




tot. sw. 


therm. 














1 


4.79 


22000 


4000 





373399 


(23) 








4.89 


20000 


2000 





386199 


(28) 








4.94 


20000 


2000 





392953 


(31) 








4.99 


20000 


2000 





400122 


(35) 








5.03 


20000 


2000 





406027 


(45) 








5.06 


20000 


2000 





410851 


(48) 








5.09 


100000 


10000 





4191 


3 


(99) 








5.0937 


530000 


30000 





4306 




(18) 








5.12 


20000 


2000 





45931 


(15) 








5.15 


20000 


2000 





46963 


(10) 








5.19 


20000 


2000 





480845 


(80) 








D.Z4 


20000 


2000 


u 


4925 


39 


(62) 








5.29 


20000 


2000 





503317 


(47) 








r - on 

5.39 


22000 


4000 





5218 


52 


(41) 








4.99 5.19 


400000 


20000 


0.402641 


(20) 


0.476307 


;29) 


0.073666 ( 


26) 


4.94 5.24 


200000 


10000 


r\ oat Ar\r\ 

0.395490 


(19) 


0.488416 


;29) 0.092926(30) 


a so k on 


200000 


10000 


n QQQ71 A 

U.ooo 1 14 


(15) 


0.499055 


V 25 ) 


0.110341 ( 


26) 


4.79 5.39 


200000 


10000 


0.375886 


(13) 


0.517555 


;2o) 


0.141669 ( 


23) 


4.99 5.03 


26000 


5000 


0.400123 


(73) 


0.405956 


'66) 






4.99 5.06 


24000 


3000 


0.400268 


(43) 


0.410796 


[72) 






4.99 5.09 


52000 


4000 


0.400280 


(37) 


0.41675 


'39) 






4.99 5.12 


24000 


3000 


0.400989 


(59) 


0.44999 


'35) 






4.99 5.15 


24000 


3000 


0.401837 


(58) 


0.46402 


;i7) 






5.03 5.19 


24000 


3000 


0.41004 


(12) 


0.47712 


;n) 






5.06 5.19 


24000 


3000 


0.41804 


(34) 


0.47823 


'19) 






5.09 5.19 


74000 


6000 


0.44597 


(67) 


0.479589 


(71) 






5.12 5.19 


24000 


3000 


0.46070 


(15) 


0.480154 


'79) 






5.15 5.19 


24000 


3000 


0.47006 


(14) 


0.48020 


;io) 







Table 2: The same as Table 1 for the 8 2 x 40 x 2 lattice. 
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ft 


run parameters 


f ) 

Pi 






Pi 




Pi- Pi 


tot. sw. 


therm. 












xlO 2 


5.6324 


12000 


2000 





533130 


(38) 






5.6524 


12000 


2000 





538 


238 


(46) 






5.6624 


12000 


2000 





540762 


(45) 






5.6724 


20000 

V 7 V 7 V 7 V 7 


2000 





543113 


(32) 






5.6< 


304 


20000 


2000 





545051 


(44) 






5.6< 


370 


20000 


2000 





54692 


(20) 






5.6924 


80000 


5000 





548 


73 


(30) 






5.6978 


20000 


2000 





552987 


(53) 






5.7044 


20000 


2000 





554618 


(51) 






5.7124 


20000 


2000 





556508 


(34) 






5.7224 


12000 


2000 





558 


488 


(36) 






5.7324 


12000 


2000 





560475 


(32) 






5.7524 


12000 


2000 





563854 


(25) 






5.6724 


5.7124 


32000 


4000 


0.544131 


(67 


) 0.555520 


'66) 


1.1390(87) 


5.6624 


5.7224 


32000 


4000 


0.541635 


(75 


) 0.557411 


; 89) 


1.5776(57) 


5.6524 


5.7324 


32000 


4000 


0.539321 


(44 


) 0.559379 


;35) 


2.0058(45) 


5.6324 


5.7524 


32000 


4000 


0.534442 


(33 


) 0.562714 


J25) 


2.8272(39) 


5.6124 


5.7724 


32000 


4000 


0.529389 


(40 


) 0.565739 


; 21) 


3.6350(47) 


5.6724 


5.682 


60000 


2000 


0.543233 


(24 


) 0.545364 


'30) 




5.6724 


5.688 


60000 


2000 


0.543296 


(33 


) 0.546843 


; 63) 




5.6724 


5.694 


100000 


4000 


0.543410 


(52 


) 0.54871 


;io) 




5.6724 


5.700 


120000 


4000 


0.543513 


(50 


) 0.55092 


; 21) 




5.6724 


5.706 


60000 


2000 


0.543827 


(41 


) 0.553736 


[77) 




5.678 


5.7124 


60000 


2000 


0.545623 


(80 


) 0.555678 


[50) 




5.684 


5.7124 


100000 


4000 


0.54831 


(11 


) o. 


556022 


[29) 




5.690 


5.7124 


88000 


4000 


0.550893 


(83 


) 0.556265 


[26) 




5.696 


5.7124 


40000 


2000 


0.552771 


(52 


) 0.556342 


[33) 




5.702 


5.7124 


40000 


2000 


0.554249 


(45 


) 0.556372 


[37) 





Table 3: The same as Table 1 for the 18 2 x 48 x 4 lattice. The last run in 
the second group was not used in the integral method. 
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N z 




5(3 




■ ojTl X 2 /dof 


0.3 


0.2 0.15 


0.1 


16 
40 


1.575(12) 
1.575 (20) 


1.084(16) 0.840(13) 
1.146(20) 0.859(20) 


0.590(14) 
0.665(15) 


0.100(18) 0.03 
0.199(34) 2.15 



Table 4: a(8(3)/T* for N t = 2 at finite 5(3 determined with the operator 
method. The last two columns are the results of a linear extrapolation to 
5(3 = and its x 2 P er degree of freedom. 



N z 




5(3 




ojTl X 2 /dof 


0.3 


0.2 0.15 


0.1 


16 
40 


1.098 (4) 
1.111(11) 


0.771 (5) 0.607(4) 
0.786(13) 0.624(12) 


0.438 (7) 
0.458(12) 


0.113(6) 0.17 
0.134(16) 0.01 



Table 5: a{5(3)/T^ for N t = 2 determined with the integral method. The last 
two columns are the results of a linear extrapolation to 5(3 = 0. 



6(3 - °ITI X 2 /dof 



0.06 om om qm 

0.6031(122) 0.3897(67) 0.2945(68) 0.2034(45) 0.0171(113) 0.06 



Table 6: cr(5(3)/T^ for N t = 4 determined with the integral method. The 
last two columns are the results of a linear extrapolation to 5(3 = using 
5(3 = 0.02 - 0.04. 
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Figure 1: The Monte Carlo simulation points for N t — 2 in the space of 
(3-). Small filled circles are the simulation points for an 8 2 x 16 x 2 lattice, 
and open circles for an 8 2 x 40 x 2 lattice. The symbol "+" corresponds to 
the phase transition point ((3 C ,(3 C ). 
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Figure 2: a with the operator method as a function of 5 ft on an 8 2 x 40 x 2 
lattice. Open circles are for our results and triangles are for those of Kajantie 
et al.\g. 
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Figure 4: Real part of the rotated Polyakov loop as a function of z on the 
8 2 x 16 x 2 and 8 2 x 40 x 2 lattices. The Polyakov loop is rotated by a Z(3) 
phase factor so that the argument is in (-7r/3,7r/3]. Errors are much smaller 
than the symbols, (a) for 5(3 = 0.1 and (b) 5(3 = 0.15. 
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d 

Figure 5: cr(d) defined as cr/T c 3 with the operator method with restricting 
the summation in (|2|) to the regions within the distance d from z = and 
N z /2. (a(N z /4) = ojTl given in Table f§) 
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Figure 6: a with the integral method as a function of 5(3. Filled circles are 
the results of the 8 2 x 16 x 2 lattice, and open circles those of the 8 2 x 40 x 2 
lattice. Triangles are the results of Huang, et al. M for N z = 16. 
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Figure 7: Plaquette expectation value in domain 1 along the path CE of 
Fig. [I] and that in domain 2 along BC on the 8 2 x 16 x 2 lattice. Estimated 
errors are much smaller than the size of symbols. 
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Figure 10: a with the integral method as a function of 6/3 on the 18 2 x 48 x 4 
lattice. The line is the result of linear fit to the data at 6/3 = 0.02 - 0.04. 
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